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Abstract 

Lattice discretization of the supersymmetric Yang-Mills quan- 
tum mechanics is discussed. First results of the quenched Monte 
Carlo simulations, for D=4 and with higher gauge groups ( 3 < 
N < 8 ), are presented. We confirm an earlier (N=2) evidence 
that the system reveals different behaviours at low and high 
temperatures separated by a narrow transition region. These 
two regimes may correspond to a black hole and elementary ex- 
citations phases conjectured in the M-theory. Dependence of 
the "transition temperature" on N is consistent with 't Hooft 
scaling and shows a smooth saturation of lattice results towards 
the large N limit. Is not yet resolved if the observed change 
between the two regimes corresponds to a genuine phase transi- 
tion or to a gentle crossover . A new, noncompact formulation of 
the lattice model is also proposed and its advantages are briefly 
discussed. 



1 Invited talk presented at the Sixth Workshop on Non-Perturbative QCD, American 
University of Paris, Paris, June, 2001. 
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1 Yang-Mills quantum mechanics and its lattice 
formulation 



The hypothesis of M-theory was a culmination of the recent developments 
in nonperturbative string theories. In the unified scheme, which emerged, 
all known five string theories are dual images of a single theory, which also 
contains eleven dimensional supergravity (for a recent review see e.g. [Q]). 
Even though it is not yet known with the full precision, the M-theory has 
a remarkable potential to unify all interactions and particles. In partic- 
ular it may offer a topological explanation of such fundamental features 
as three families and fractional charges. It may lead to a standard model 
gauge group with J\f=l supersymmetry. It provides understanding of the 
Bekenstein-Hawking entropy puzzle and much more. Due to the proposition 
of Banks, Fishier, Susskind and Shenker j|| , M-theory became amenable to 
yet more quantitative studies. According to BFSS the spectrum of M-theory 
is equivalent to that of a super symmetric Yang-Mills quantum mechanics 
(SYMQM). The latter is obtained from the dimensional reduction of the 
10 dimensional supersymmetric Yang-Mills theory to the one (time) dimen- 
sion. Such a quantum mechanical system can be quantitatively studied by 
a host of nonperturbative methods. One of them, which provides steadily 
increasing understanding of nonperturbative phenomena, is the lattice ap- 
proach. Therefore we have constructed the Wilson discretization of the 
above quantum mechanics and proposed to investigate it with the standard 
lattice methods ||. The ultimate goal is to study D=10 SYMQM for the 
large size of the SU(N) matrices. However, even for this relatively simple, 
one dimensional quantum mechanical system this is still a rather complex 
task. One of the main difficulty is posed by the complex fermionic determi- 
nant (and pfaffian) for addjoint fermions in D=10. Another, although less 
severe one, is the time consuming Monte Carlo calculation with large non- 
abelian matrices. On the other hand the system can be simulated at present 
in other regions of the parameter space D, N, Nf, (Nf being the number of 
fermions ) and such a study may provide a relevant information about its 
general structure. Therefore we have decided to set up a systematic lattice 
survey of SYMQM beginning with the simplest case of D = 4, N = 2, Nf = 
and gradually extending it as far as possible towards the BFSS limit. In this 
talk I will report on the second step along this program, namely the first 
results for higher N will be presented. 

Supersymmetric Yang-Mills quantum mechanics [Q, ||, |(| and its zero 
dimensional counterpart 0, || have been intensively studied. Although the 
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exact solution is still not available, many results are known and a) can be 
tested and extended, and b) provide us a guidance from a "simple corner" 
of parameter space to the ultimate "BFSS corner". 
The action of the SYMQM reads 

S = J dt QttF^) 2 + ^{t)V^ a {t)\ . (1) 

where fi, v = 1 . . . D, and all fields are independent of the space coordi- 
nates x. The super symmetric fermionic partners belong to the addjoint 
representation of SU(N). The discretized system is put on a D dimensional 
hypercubic lattice N± X . . . X Np which is reduced in all space directions 
to Ni = 1, i = 1 . . . D — 1. Gauge and fermionic variables are assigned to 
links and sites of the new elongated lattice in the standard manner. This 
geometry has two straightforward consequences. First, all space derivatives 
in bosonic and in fermionc terms vanish, di — > 0. Second, there is a new 
class of gauge invariant observables, namely Tr(Ai{t) 2 ) or Tr(C/j(m)), since 
a local gauge transformation of the spatial components A; L (t) reduces to a 
simple similarity transformation. The gauge part of the action has the usual 
form 

S G = -PJ2J2 ^e(TrtVM), (2) 

771=1 ii>V 

with 

P = 2iV/aV, (3) 

and U^yim) = Ul(m)U^(m+v)U v (m+ii)U^(m), U^m) = exp (iagA^am)), 
where a denotes the lattice constant and g is the gauge coupling in one 
dimension. The integer time coordinate along the lattice is m. Periodic 
boundary conditions U^{m + v) = U^(m), u = 1 . . .D — 1, guarantee that 
Wilson plaquettes tend, in the classical continuum limit, to the appro- 
priate components Fa V with the space derivatives absent. 

Another simple property of the one dimensional gauge system shows up 
in Eq.(||). The gauge coupling g has a dimension, and as a consequence 
its numerical value provides directly a scale to all phenomena occuring in 
this system. In a language familiar from the four dimensional QCD: the 
dimensionless lattice coupling (3 = j3{a) is running according to Eq.(||) , and 
g plays a role of A QC d- 

In the above formulation the projection on the gauge invariant states is 
naturally implemented. 
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Eq. (^) was the basis of the MC simulations for a simplified SU(2) model 
(see Sect. (2.1)). Now some results for higher gauge groups N < 9 are also 
available and will be discussed in Sect. (2.2). 



2 Results 
2.1 SU(2) 

One of the remarkable features of the M-theory is the explanation of the 
problem of the Bekenstein-Hawking entropy || in terms of the elementary 
excitations. Moreover, the theory predicts existence of at least two phases: 
a low temperature, "black hole", phase and a high temperature phase de- 
scribed by the elementary DO branes 0. Therefore a natural question 
emerges if the simplified model (i.e. SYMQM), considered at finite temper- 
ature, possesses any nontrivial phase structure [11]. It is well established 



that QCD (or pure Yang-Mills theory) has two different phases (e.g. con- 
finement and deconfinement). Since the action (||) is basically QCD-like, 
one might naively expect that the dimensionally reduced model may indeed 
exhibit similar phenomenon. On the other hand the one dimensional model 
with local interactions cannot have a phase transition for finite N. However 
at infinite N the sharp singularity may occur [O]. It follows that, contrary 
to the first expectations, the mechanism which may eventually lead to such 
a transition will have to be rather different that that for, e.g., N = 3 QCD. 
Hence the problem is even more interesting, since a) it is relevant for the 
M-theory builders and b) the nature of the transition (if any) is open. 

In Ref.[|| we have studied the phase structure of the quenched model 
with D=4, N=2 and Nf = 0. Since N plays a role of the volume, we 
expected to see some signatures of a phase change for finite and even small 

To check this we have measured the distribution of the trace of the 
Polyakov line 

P = ^ Tr (li PdH) • ( 4 ) 

which is a very sensitive determinant of the two phases in QCD. Similarly to 
lattice pure gauge system, symmetric concentration of the trace around zero 



2 The full phase structure of the M-theory is much more complex, op. cit. 

3 It is well known from statistical physics that the singularities (in the temperature) 
associated with the phase transition, develop gradually and smoothly with increasing the 
volume and show up ( as a broad peaks for example) even for small systems [p.4|. 
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indicates the low temperature phase with < P >= (here a "black hole" 
phase) while clustering around ±1 (or, for arbitrary N, around the elements 
of Zn) is characteristic of the high temperature phase (here the elementary 
0-brane phase). It follows from Eq.© and N t = l/Ta that = 2NN?T 3 , 
hence the lattice coupling provides, at fixed N t , a direct measure of the 
physical temperature T . It was found that the distribution of (|j) is indeed 
changing from a convex to a concave one at some finite value of (3 = (3 c {Nt) ■ 
Moreover the dependence of (3 C on Nt is consistent with the canonical scaling. 
The best fit to Monte Carlo results gave (3 c {N t ) = (0.17 ± 0.05)ivf - 02±0 - 33 ) 
in good agreement with the expected dependence (3 C ~ A^ 3 . This means that 
in the continuum limit the transition occurs at finite temperature T c , which 
can be written as T c = ^(c^A^ 1 / 3 ) with A 2 = .28 ± 0.03. 

We have also measured the average size of the system R 2 = g 2 J2 a ( A i f 
for different values of the temperature. It shows entirely different behaviours 
at low and high temperatures with a narrow transition region. The pseu- 
docritical temperature determined in this way is consistent with the one 
obtained above from the study of the Polyakov line. Our results also agree 
qualitatively with the mean field calculation at large ./V [0|. Obviously 
there is a long way between N = 2 and the BFSS limit N — > oo and it was 
essential to extend this study to the larger groups. 



2.2 SU( 3 < N < 8) 

Recently we have done Monte Carlo simulations with the intermediate size 
groups, still for Nj = and D = 4 [jDJ. We used the standard Metropo- 
lis algorithm employing all independent moves in the SU(N) manifold. For 
N > 2 the trace P in Eq.(|j) is complex. The condition det(U)=l constraints 
P to lie inside the "N-star" in the complex plane (c./. Fig.l). New simu- 
lations entirely confirm the structure found in the SU(2) case. In the low 
temperature phase the distribution is symmetric and peaked around P = 
which results in < P >= 0. At high temperatures P's are expected to con- 
centrate around the elements of Z^. A sample of our results for N=3,5,8 
and at high and low temperatures (/? ~ T 3 ) is shown in Fig.l. Indeed 
one clearly sees the existence of both regimes. Moreover, in the high tem- 
perature region, the system has a tendency to be stuck in one of the Zjy 
minima. Correspondingly the thermalization time grows strongly with N. 
This is a typical indication of the spontaneous symmetry breaking which 
might occur at N — > oo leading eventually to a non-zero value of the order 

4 Up to the finite a corrections. 
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Figure 1: Distribution of the Polyakov line in the low and high tempera- 
ture phases (upper and lower rows respectively) of the quenched Yang- Mills 
quantum mechanics with the SU(3), SU(5) and SU(8) gauge groups. 

parameter < P >. 

We cannot determine at present the nature of the transition separating 
the two phases. In fact it is also possible that this is not the phase transi- 
tion at all, but just a gentle crossover between the high and low temperature 
regimes. More extensive simulations are required to answer this question. 
Nevertheless, our data show unambiguously that the system behaves differ- 
ently in both regions. For example, the dependence of the size of the system 
on the temperature is definitely different, and the change in the behaviour 
occurs at the same T where the distributions in Fig. 1 change their shapes. 
Moreover, the characteristic temperature separating the two regions (loosely 
called phases) is physical in the sense that it is finite in the continuum limit. 
This follows simply from the scaling of the lattice coupling f3 c . 

With present data one can study, for the first time, the N dependence of 
the transition temperature. It follows from Eq.(||) and Nt = 1/Ta, that the 
't Hooft scaling T c ~ (g 2 N) 1 ^ 3 implies that the lattice coupling /3 C ~ N 2 at 
fixed Nf. This prediction is tested in Fig. 2, where (3 C /N 2 is plotted versus 
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Figure 2: Dependence of the transition lattice coupling on the number of 
colours N at fixed time extent iVj = 4. 



N. The values of [3 C were determined for each N, analogously to Sect. 2.1, 
as the couplings for which the distributions, Fig.l, change their shape. The 
errors represent the subjective uncertainty of such an estimate. Indeed a 
nice saturation of the N dependence is observed in Fig. 2 confirming the 't 
Hooft scaling of T c . The solid line is a fit of the 1/N 2 correction, and the 
flat solid line indicates the fitted asymptotic value. The data for the first 
two N require yet higher order corrections, while simulations for SU(8) give 
the asymptotic value with less than 15% accuracy. 

It is interesting to compare above results with the N dependence found in 
the full, spatially extended Yang- Mills theory [15] where the 1/N 2 correction 
seems to work well even for N=2,3 in some cases. It would be instructive to 
repeat their study for the finite temperature phenomena. 



3 Noncompact formulation 

In this section we propose an alternative discretization of the continuum 
model (|l]) , which turns out to be easier to simulate numerically. The idea is 
to leave the D — 1 spatial components of the gauge fields noncompact, while 
keeping the D-th component compact. To this end split the continuous 



action of the reduced theory into the kinetic and potential terms. 

S = S^in + Spot, 

S km = jzjdtTriDtX') 2 . 



where the sums over i and k are implied and D t denotes the covariant 
derivative along the time direction After discretizing the time we obtain 

Spot = 2?Em^ r (^ , ^ t ) 2 ' 

Skin = ^E m Tr(AX^, (6) 

with the D — 1 noncompact matrix coordinates X*(m) defined at a discrete 
time intervals t m = am. The covariant finite difference along the time 
direction 

AX'(m + 1) = X*(m + 1) - U(m + 1, m)X i (m)C/(m, m + 1). (7) 

takes into account the parallel transport between adjacent lattice cites. The 
system now has D — 1 noncompact coordinates X* and one compact degree 
of freedom corresponding to the timelike link U(m + l,m) = £/£>(m). The 
discretized action is invariant under the same local gauge transformations 
as the compact version (^) 

X l (m) -> V- 1 (m)X. t (m)V(m), 
U(m + l,m) -> V- 1 (m + l)U(m + l,m)V(m), (8) 



As a consequence of the dimensional reduction there is no inhomogeneous 
term in the gauge transformation of the spatial coordinates X\ Also, they 
are rotated by the same element of the gauge group V(m) on both sides, since 
the space like links on Ni = 1 lattice , Eq.(^) close back onto themselves. 

With the noncompact action (JgJ) one can readily perform simulations 
closer to the continuum limit than those done with the fully compact vari- 
ables. Also higher N's can be easier achieved. This is because the lattice 
coupling = 2N/a s g 2 diverges strongly at small a and this causes very 
severe critical slowing down in the compact case. The problem is even more 
difficult for larger N. On the other hand the effective coupling in the non- 
compact formulation (|6|) diverges only as \ja which alleviates the above 
difficulty. 
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We have done some exploratory simulations with the new formulation for 
a range of iV's (2 < N < 9) and for different dimensions of unreduced theory 
(4 < D < 10). All simulations in all these dimensions confirm existence of 
the two different regimes of the model in full analogy with the compact, 
N = 2 case. We have encountered no difficulty extending these study to 
higher D. In fact we have found a steady decrease of the average size of a 
system (defined per one link) with D, in both phases, in agreement with the 
mean field predictions 

4 The future 

Lattice simulations provide a new approach to a quantitative study of the 
Yang-Mills quantum mechanics and possibly to the M-theory. Preliminary 
results are encouraging, but a lot remains to be done. Simulations of the 
quenched model work for all interesting values of the dimension D and are 
feasible for a range of N. Recent results give us a rough idea how the large 
N limit is approached and where the asymptotics sets it. All quenched 
simulations performed up do date indicate existence of the two regions at 
finite temperature. This intriguing correspondence with the predictions of 
the M-theory should be (and can be) further quantified. Of course, the next 
step is to include the dynamical fermions. This can be done by a brute 
force for D=4 and for first few N's at D=10. The one dimensional nature 
of the system should help considerably. For higher N, at D=10, we face the 
problem of the complex pfaffian and new ideas are needed. Finally, the full 
potential of the small volume approach [fL6j |l7| may provide an important 
insight into the problem and should be applied to these supersymmetric 
systems. 
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